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Abstract. We here report calculations of magnetic buoyancy instabilities of a sheared magnetic layer 
where two separate zones are unstable. The idea is to study the possible generation of large-scale helical 
structures which could then rise through a stellar convection zone and emerge at the surface to create active 
regions. The calculations shown here are a follow-up of the work ofjFavier et al.| ( [2012 ) where the instability 
developed in a weakly magnetized atmosphere, consisting of a uniform field oriented in a different direction 
from the unstable layer below. Here, the top layer representing the atmosphere is itself unstable to buoyancy 
instabilities and thus quickly creates a more complex magnetic configuration with which the layer below will 
interact. We also find in this case that the accumulation of magnetic tension between the two unstable layers 
favors the creation of large-scale helical structures. 
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1 Introduction 

The idea of this work is to study the onset and development of magnetic buoyancy instabilities of a sheet of 
field in a configuration where large scale helical magnetic structures may emerge. Detailed calculations were 



performed by Favier et al. (2012), in which the buoyancy instability develops in a weakly magnetized atmosphere 
where the direction of the field varies with respect to the direction of the strong unstable field below. The idea 
was to understand how coherent structures could be created at the base of a stellar convection zone where a 
large-scale shear produces a strong toroidal magnetic field associated with a small poloidal component. The 
authors showed that from small-scale interchange instability, large-scale helical coherent magnetic structures 
could be created, thus giving a possible explanation for the coherent rise of magnetic flux from the base of the 
convection zone to the surface where bipolar active regions emerge. 

Here, we investigate the case where the atmosphere above is itself unstable to buoyancy instabilities and 
naturally creates a more complex field configuration. The idea here is to study the possibility of the spontaneous 
creation of similar coherent structures when the atmospheric field is less idealized and organized at smaller 
scales. The calculation thus starts with an initial magnetic layer oriented in the x-direction, complemented with 
a transverse field By such that two zones are initially made unstable to 2D and 3D perturbations. 

2 Model set up 

2.1 Governing equations 

The evolution is governed by the compressible MHD equations, with heat transfer treated in the diffusion ap- 
proximation. The polytropic atmosphere is considered to be convectively stable so as to focus on the instabilities 
driven by the magnetic field. Those equations are solved in a Cartesian domain in the same way as |Matthews] 



et al. (1995). In particular, the following dimensionless parameters are used: 
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where po^ po and Tq being the pressure, density and temperature at 2: = 0, /io the magnetic permeabihty, Cp 
the specific heat at constant pressure, R the gas constant, g the gravitational acceleration, A is the temperature 
gradient and d is the depth of the box. The thermal conductivity viscosity u and magnetic diffusivity r] are 
all assumed to be constant. 



The equations then read: 



g = -V-(pu) (2.1) 



^ = -V(p + FB^/2) + V ■ (FBB - puu + aC'kT) + (m + l)0pz (2.2) 

^ = Vx(u X B - CCfcV X B) (2.3) 

= -u • VT - (7 - 1)TV • u + l^V^T + ^'^^~^^ (arV2 + FCJ') (2.4) 
at p p 

where r is the stress tensor and where the ideal gas law is used. The boundary conditions used are imper- 
meable and stress-free at top and bottom, vertical magnetic field at top and bottom and constant fiux at the 
bottom for the temperature. 

2.2 Initial field configuration 

The initial magnetic field configuration consists in a layer of magnetic field in which two spatially separated 
zones will be unstable with respect to buoyancy instabilities. The magnetic layer is located between z = 0.6 



and z = 0.8 (the top of the domain being located at z=0). As shown by Newcomb (1961), the criteria for 2D 
and 3D instabilities respectively read: 



dz \ p J dz 

and 



The profile of the transverse field is chosen here so as to make the left-hand sides of Eq 2.5 and 2.6 positive 
in two different locations between z = 0.6 and z = 0.8. As a consequence, both the overlying atmospheric field 
and the layer below (possessing different field inclinations) will be unstable to buoyancy instabilities. 

As initial profiles for the magnetic field we choose the following: 

= Bo{z- 0.6)"^ {z - 0.8) 

By = Boy {z - 0.6)^/^ {z - 0.8) exp 

Figure [l] shows the initial magnetic field profile as well as the characteristic quantities for the onset of magnetic 
buoyancy instabilities. It is clearly visible on the right panel of this figure that the magnetic field profile is 
made unstable at 2 different locations within the layer, without any interaction between them at the initial 
state. The profiles of B^ and By are smooth, the longitudinal field keeps the same sign within the whole layer 
while By changes sign at z = 0.7. Of course, in this case, the top unstable layer is dominated by a field in the 
^/-direction, which would correspond to an almost entirely poloidal field at the base of the convection zone, and 
thus probably unlikely. Nevertheless, we are interested in the effect of a transverse field varying with depth and 
this configuration will thus capture the largest possible effects of this situation. 

The dimensionless numbers used in this simulation are the following: 7 = 5/3 , m = 1.6,0 = 2, C/c = 
1.25 X 10~^, cr = 2 X 10~^ for the hydro numbers and F = 1, = 2 x 10~^ for the magnetic quantities. The 
resolution used is = 256, Ny = 128 and = 260. 
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Fig. 1. Left panel: initial profiles of Bx and By. Right panel: initial values of the right hand side and left-hand sides of 
Eqs [2:51 and 



3 Results 

3.1 Layer interaction 





Fig. 2. Evolution of By in the two unstable zones within the domain. By is positive in the top layer (left panels) and 
negative in the bottom layer (right panels). The isosurfaces have dimensionless values between —0.6 and 0.6, with a 
distance of 0.1 between them. The colour scale is thus kept the same for all snapshots, allowing us to see the effect of 
magnetic diffusion on the rising layers. 

Figure |2] shows isosurfaces of the two unstable layers in this case. We focus on each unstable zone separately, 
in order to have a clearer idea of what the orientation of the field is in each layer and how the instability sets 
in. The two different orientations of the field appear clearly in this picture. We note that the positive By 
zone (the top unstable layer) gets unstable before the other layer and very narrow flux tubes start to emerge. 



This is consistent with the linear analysis of Cattaneo et al. (1990) and Favier et al. (2012) which showed that 



the presence of an overlying atmosphere with inclined field lines reduces the growth rate of the instability of a 
magnetic layer. 

The two layers then interact and create 3D structures. At t = 14.89, the structure of the top unstable layer 
is very interesting, small arched flux tubes are created and tend to maintain good coherence, a configuration 
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which is quite favorable for the emergence of twisted structures. This is indeed what happens, as wih be shown 
in the next subsection. Both layers now emerge with different speeds and rise towards the top of the layer. We 
find here that the top unstable layer rises on average more slowly than the bottom one. It is indeed affected 
more quickly both by viscosity and magnetic diffusion since the instability grows faster and since the flux ropes 
are smaller scale. As a consequence, at around 20-30 sound crossing times, the unstable modes of the bottom 
layer start to interact with those of the upper layer and the imprint of the preferred orientation of the modes 
of the blue layer can be seen on the red layer (top one). There is thus a strong interaction between the modes 
of various orientations and this will have consequences on the field line structure. 

It is worth noticing that the whole layer gets completely disrupted by the magnetic buoyancy instabilities in 
this case, even though the resonant layer is located in the middle of the magnetic domain (see discussion on the 



resonant layer by Cattaneo et al. 1990). The vertical velocity (not shown here) indeed shows that downward 



motions are very strong, concentrated on narrow bands and penetrate very deep into the magnetic domain. It 
thus results in the complete disruption of the magnetized zone after a few tens of sound crossing times. 




Fig. 3. Cuts at 2/ = est of the magnetic field components Bx and By. For Bx, the colour scale ranges between and 0.8 
and for By, the range is —0.6/0.6. Large-scale magnetic structures are clearly visible on both snapshots. 



Figure [3] shows cuts at y = est and at two specific times of the x and y components of the magnetic field. 
This picture shows the interaction between the very small scales emanating from the interchange instability 



of the top layer and the much larger scales below. As shown by Favier et al. (2012) in their more idealized 



atmosphere, the magnetic tension lying at the interface between the two unstable layers promotes the merging 
between small-scale structures in the bottom magnetic layer. It is indeed quite clear that between t = 11.17 
and t = 22.32, the top layer has been completely disrupted into small-scale features when the bottom layer has 
built much larger scales. 

As a consequence, large-scale coherent structures are also generated in this less idealized case and could in 
principle rise through a large number of pressure scale heights before being diffused away or reprocessed by 
convective flows in a real stellar convection zone. 
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3.2 Generation of twisted flux ropes 



Finally, since our initial goal was to study the possible generation of twisted flux ropes in these calculations, 
we focus on the field line configurations, especially on the small unstable structures emerging in this simulation 
at the top of our magnetic domain. Figure [4] shows some rendering of magnetic field lines around current 
helicity concentrations in this calculation. As also shown by Favier et al. (2012), the initial configuration 



chosen here not only favors the generation of large-scale magnetic features from small-scale magnetic buoyancy 
instabilities but also promotes the creation of twisted flux ropes. Those helical fiuxtubes should evolve in the 
polytropic atmosphere keeping much more coherence than the magnetic structures which are produced in a 
purely horizontal field case. Indeed, it was shown in multi-dimensional simulations of the dynamics of magnetic 
flux ropes (e.g. Emonet fc Moreno-Insertis]|1998 ) that a certain amount of twist of the field lines was necessary 
for concentrations of magnetic field to compete against the gravitational torque exerted on the rope during its 
emergence. Figure [4] shows that in these calculations, helical magnetic structures naturally emerge from the 
buoyancy instability in an initial configuration where both a longitudinal and a transverse field are present. 




Fig. 4. Field line tracing around small flux ropes created in this simulation through the instability of a magnetically 
sheared layer. 



4 Conclusion 



In this work, we studied the development of magnetic buoyancy instabilities of a magnetic layer in which two 
zones with different orientations of the field lines become unstable. The top layer becomes unstable first, in 
agreement with linear theory, and the two layers quickly interact, creating large-scale helical structures at the 
interface where magnetic tension and twist accumulate. The outcome of those calculations is similar to the case 
studied by [Favier et al. (2012) where a layer of field becomes unstable and interacts with a weakly magnetized 
(and thus stable) atmosphere. We thus show here that despite the more complex configuration of the atmosphere 
above (created here by the unstable top layer), twisted large-scale magnetic structures can still naturally emerge 
from buoyancy instabilities. This could in principle explain how strong concentrations of field created at the 
base of the solar (and stellar) convective zones are able to rise coherently to the surface to create the observed 
active regions. 
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